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Abstract 

It is shown that a strong intraplanar incoherent scattering can effectively 
block the interplanar coherent tunneling between the weakly coupled planes 
of the highly anisotropic but clean (intrinsic) materials such as the op- 
timally doped high-Tc layered cuprate superconductors. The calculated 
normal-state C-axis resistivity Pc{T) then follows the metal-like tempera- 
ture dependence of the ab-plane resistivity PabiT) at high temperatures. At 
low enough temperatures, however, Pc{T) exhibits a non-metal like upturn 
even as Pab{T) remains metallic. Moreover, in the metallic regime, Pc{T) is 
not limited by the maximum metallic resistivity of Mott-Ioffe-Regel. This 
correlation between the intrinsic Pc{T) and Pab{T) is observed in the normal 
state of the high-Tc stoichiometric cuprates. 
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1. Introduction. 

The normal-state out-of-plane resistivity Pc{T) of the high-Tc layered 
cuprate superconductors^ has raised a number of questions, ^'^ as yet un- 
resolved, about its temperature dependence, its absolute magnitude, and 
its dependence on the concentration of carriers, i.e., doping. Thus, while 
the in-plane resistivity Pab{T) is well known to be metallic^ (i.e., with 
TCR = dpab/dT > 0, and, in fact, essentially T-lincar right from Tc up- 
wards to the highest temperature of measurement) and smaller than the 
Mott-Ioffe- Kegel maximum metallic resistivity p^^x^ out-of-plane re- 
sistivity Pc{T) shows a range of behaviour. Thus, it has been variously 
reported to be non-metallic^~^ {TCR < 0) for underdoped samples; mixed- 
met allic^'^'^^ (metallic at high temperatures but with a non-metallic uptrun 
at low enough temperatures); and completely metallic^"^'^"^^ for stoichio- 
metric (fully oxygenated) composition showing a T- linear Pc{T) from Tc 
upwards, but with the absolute magnitude of pdT) > p^ax cases. 
The essential structural feature of weakly coupled layers and the associated 
large resistive anisotropy with pc{T)/pab{T) ~ 10^ — 10^, clearly makes the 
C-axis resistivity highly sensitive to extrinsic details that presumably tend 
to contaminate its intrinsic behaviour. Thus, it is entirely possible for any 
measurement of the out-of-plane resistivity to pick up some in-plane com- 
ponent of the resistivity tensor — perhaps externally due to misalignment 
of the contacts, or internally due to the randomly distributed defects and 
faults providing shorts between the weakly coupled ab- planes. Such a con- 
tamination of the out-of-plane pdT) by the in- plane Pab{T) can make Pc(T) 
track the metallic temperature dependence of Pab{T), making the former ( 
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Pc{T)) an apparent metal}^ The resistivity data on high quahty untwinned 
single crystals, however, strongly suggests that the out-of-plane resistivity 
Pc{T) is intrinsically metallic, and in fact T-linear^^'^^ at least at high tem- 
peratures, just as the in- plane pabiT) is. Its absolute magnitude is, however, 
much larger. Recently, we had proposed a mechanism that gave precisely 
such a behaviour. In this mechanism^^ the inter-planar tunneling between 
the weakly coupled metallic planes is cut-off (blocked) by the intra-planar 
inelastic (incoherent) scattering, leading to PciT) oc PabiT). This physical 
mechanism has since gained a fair degree of acceptance among the workers 
in the field, ^ while some earlier theories, linked closely to the exotic mecha- 
nisms for high-Tc superconductivity in the strongly correlated Cu02 sheets, 
have been argued out to be inconsistent with known experimental facts. ^ 
Motivated by these developments, we have re-examined the mechanism pro- 
posed by us earlier based on simple physical arguments. In doing so we have 
derived an expression for pc{T) following the Kubo-Matsubara conductivity 
formalism applied to a model Hamiltonian incorporating weak interplanar 
tunneling and strong intraplanar incoherent scattering. Our Pc{T) so de- 
rived indeed shows a metallic behaviour with pc{T) oc Pab{T) in the high 
temperature limit, thus validating our mechanism proposed earlier. At low 
enough temperatures, however, we get an additional feature of a resistivity 
upturn {dpc/dT < 0), which is qualitatively consistent with observations as 
noted above. In the following, we give some details of our derivation, and 
discuss our results for Pc{T) in the light of some recent findings. 

2. Theoretical. 

First note that the C-axis transport, except possibly for the overdoped 
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cuprates, is known to be incoherent, e.g., ojpcTab <S 1, where ojpc is the 
C-axis plasma frequency and 1/Tab is the intraplanar inelastic scattering 
rate.2'^^ Thus, the successive interplanar tunneling amplitudes are phase- 
uncorrelated. It is, therefore, sufficient to consider simply a bilayer (AB) 
coupled weakly by a tunneling matrix element (— tc)- In real systems the 
individual layer (A or B) can by itself represent a single Cu02 sheet as 
in LSCO, or also a group of strongly coupled Cu02 sheets, as in YBCO, 
BSCCO and other multilayered cuprates, separated by the spacer oxide 
layers. Also, we will consider only the clean hmit as suggested by the 
smallness of the zero-temperature intercept^^ PabiT — > 0), and also assume 
that the inter-planar tunneling conserves the wavevector parallel to the ab- 
plane. Then the model Hamiltonian (in obvious notation) is 

H — Ha + Hi, + Hab + Haa + -f^&6 , (1) 

with 

Ha - y: 

kcr 

Hb = Xl^k&L^ka 

kc 

Hab = -^cX!(4a^ka + &L«ka) 
k(7 

In the following we will drop the spin index cr. Here Haa and i^bb represent 
the inelastic intra-planar electron- electron scattering characteristic of the 
strongly correlated two-dimensional Cu02 sheets. In the present analysis, 
however, these terms shall enter only implicitly and summarily through the 
imaginary part of the associated retarded electron self-energy chosen so as 
to be consistent with the known T-linear in-plane resistivity PabiT)}^ Our 
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problem then is to calculate the out-of-plane resistivity Pc{T), given the 
above PabiT) as an input. It is really this connection between PciT) and 
Pab{T) that is being addressed here. 

Now, for the inter-planar current operator, we have 

jc = -ietc (4^k - blttk) . (2) 



The Kubo conductivity in the dc hmit is then given by {h — 1) 
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a = -(-V-0 , (3) 

where c = bilayer separation (the C-axis lattice constant) , — area of 
the layer, and the retarded correlation 

Uretiuj) = limn(icu^) (4) 
icu,^ cu + iS (analytic continuation) 

with Il{iu!^), the current-current correlation given by 

1 

= 2eHl ^ -5Ga(k, ia;„)Gb(k, iuj^ + ioJv) (5) 
k 

with 

ujv = the Bosonic Matsubara frequency, 

ojn — ^ — ^— , the Fermionic Matsubara frequency. 

/-^ 

Here Ga and Gb are the temperature Green functions for the layers A and B, 
respectively, in the presence of inter-planar tunneling. For identical layers, 
as in the present case, we have Ga = G^, = G, say. 

In the high temperature limit, i.e., for h/Tab ^| tc |, we can evaluate G 
in the presence of interplanar tunneling in terms of Go, the corresponding 



temperature Green function in the absence of tunneling, from the Dyson 
equation 

Ga — Goa + Goa tc Gbtc Goa ■, (6) 

giving 

G 

Ga ^ Gb = G ^ ^ __^2 Q2 ■ ^'^^ 

In writing Eqn. (6) we have used the approximation that for h/rab ^| tc | 
one can neglect the vertex correction arising from the interplanar tunneling 

tc- 

Now, the intra- planar thermal Green function Go(k, iu!n) corresponds to 
an isolated layer (i.e., with tc — 0) and has the general form 

lUn - ek - E(k, lUJn) 

Substituting from Eqns. (8) and (7) into Eq. (5), we get 

Tl{iujy) ^ 2eHl ^ G„^(k,ia;n)Go^(k,ia;„ + ia;,,) (9) 

k:,n 

Here ± refers to Go in Eq. (8) with ek replaced by ek± \tc\. 

Now we impose our condition of T-linearity of Pab{T) as input at the level 
of Go(k, iuJn) namely, that the self-energy of the corresponding retarded 
Green function must have an imaginary part A(T) (at the Fermi level) 
= —ImJ2ret ^ T. With this input Eq. (9) together with Eq. (3) gives, 
after the usual frequency summation and analytic continuation, a simple 
expression for the dc conductivity 

Here we have introduced the two-dimensional density of states assumed 
constant. Also, fi has been re-instated. This is our main result. lYmm. 



Eq. (10) it is readily seen that in the high-temperature Hmit, A(T) 3>| tc I 
, cr oc 1/(A(T), or equivalently, the C-axis resistivity Pc{T) oc A(T) oc T, 
confirming the T-hnearity of Pc{T) at high temperatures. It is also clear that 
this mechanism giving incoherent transport along the c-axis does not involve 
the usual '/ci^^'-parameter characteristic of metallic transport. Hence, Pc{T) 
is not subject to p^ax- 

Next, wc consider the low-temperature limit, h/Tc <\ tc |. Now, we must 
diagonalize the tunneling Hamiltonian Ht = Ha + Hf, + Hab first, and then 
use the T-linearity of Pab{T) as an input at the level of the layer-diagonal 
Green function. We have 

Ht = Z!(^k- I tc \)al ftk + Z!(^k+ I tc \)/3l /3k , 

k k 

where = (ok + &k)/\/2, Pk = (ok — &k)/\/2- The inter-planar current 
operator 

j = ietc J2 i^l /^k - PI oik) ■ 

k 

Repeating the earlier steps with this current operator, we now get 



Thus, at low enough temperatures we get an upturn for pc because A(T') oc 
T. This upturn has been noticed as discussed earlier. It must be emphasized 
here that this upturn is a consequence of our assumption of conservation 
of wavevector parallel to the planes in the tunneling process. This leads 
to a hybridization gap that suppresses the overlap of the spectral functions 
corresponding to Go+ and Go-- 
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3. Discussion. 

We would now like to comment on a recent attempt to explain away 
the metallicitiy of Pc{T) as a thermal expansion effect. ^° Based on a recent 
study of the pressure dependence of Pc{T) and the lattice constant 'c' in 
La2-xSrxCu04, it has been argued that the observed metallicity (i.e., TCR 
> 0) merely reflects a thermal expansion effect, i.e., 

TCR = idpJdT)p = {dpJdP)T idc/dP)^^ {dc/dT)p , 

and, therefore, the observed metallicity is only apparent. While some such 
expansion effect cannot be ruled out, we find the argument somewhat flawed 
for the following reason. The thermal expansion < 5c > of the lattice pa- 
rameter 'c', or more precisely < SB > of the width 'B' of the potential bar- 
rier for the rate-determining tunneling, is a mean-anharmonic effect that, of 
course, tends to diminish the tunneling rate. There is, however, also a ther- 
mal fluctuation about this mean- value which is of comparable magnitude, 
if not larger. Now, inasmuch as the tunneling time is expected to be much 
shorter that the typical lattice vibrational time period, we must average the 
instantaneous tunneling rate over the anharmonic fluctuations of SB. As 
the tunneling rate depends exponentially on SB, overall it is expected to 
produce an enhancement of the tunneling rate with increasing temperature, 
and hence a negative, rather than positive TCR. 

Finally, a remark on the recent attempts^^'^^ to explain the resistivity 
upturn in terms of pre-existing real-space pairs (Bosonic) as precursor to 
their condensation at Tf.. Such a system of charged Bosons in the normal 
state is, however, expected to exhibit large, universal diamagnetism which 
is not reported. 
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In conclusion, the present microscopic-treatment of the C-axis resistivity 
supports the mechanism proposed by us earher, namely that the strong 
intra-planar incoherent scattering cuts-off the interplanar tunneling, and 
thus correlates Pc{T) with pabiT). 
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